Mesoscopic fluctuations of nonlinear conductance of chaotic quantum dots 
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The nonlinear dc conductance of a two-terminal chaotic cavity is investigated. The fluctuations 
of the conductance (anti)symmetric with respect to magnetic flux inversion through multichannel 
cavities are found analytically for arbitrary temperature, magnetic field, and interaction strength. 
For few-channel dots the effect of dephasing is investigated numerically. A comparison with recent 
experimental data is provided. 

PACS numbers: 73.23.-b, 05.45.Mt, 73.21. La, 73.50.Fq 



X 



Introduction. Recently the non-linear dc conductance 
of mesoscopic structures and, more specifically, their 
dependence on magnetic flux $ has found considerable 
attention 0, 0, 0, El ■ Application of large voltages in- 
duces a rearrangement of the charge distribution. The 
charge re-distribution is subject to Coulomb interactions 
0. Consequently investigation of non- linear transport 
reveals information on interaction parameters. This is 
in marked contrast to linear transport where the dc con- 
ductance G a p = dla/dVp for large cavities can with high 
accuracy be treated within a theory of non-interacting 
electrons. However to extract this information the role 
of temperature and dephasing need to be known with 
precision. 

Interaction constants are extracted by investigating 
the magnetic field symmetry of non-linear transport. Un- 
der flux <!> reversal, in the linear regime, the Onsager- 
Casimir relations dictate that the conductance matrix 
has the symmetry G a p{<&) = Gp a {— In particular, 
the conductance of a two probe conductor is an even func- 
tion of flux Q. However, away from equilibrium, the non- 
linear conductance lacks such a symmetry. Importantly, 
the deviations from Onsager symmetry are entirely due 
to interactions 0, Ej • Therefore by investigating the de- 
parture from the Onsager-Casimir relations, information 
on the interaction properties can be obtained. 

Our work is motivated by very recent experiments on 
nonlinear transport in various open systems: carbon nan- 
otubes 0], quantum dots 0, ballistic billiards Q, and 
quantum rings ||. Because of quantum interference, the 
samples exhibit strong mesoscopic (sample-to-samplc) 
fluctuations, and a theory has thus to predict statisti- 
cal properties. Only recently have two theories explored 
such statistics for two-terminal open samples: Sanchez 
and Buttiker |8( considered chaotic quantum dots in the 
universal regime with arbitrary interaction strength and 
high magnetic fields at T = 0, and Spivak and Zyuzin 
concentrated on weakly interacting open diffusive sam- 
ples at low fields and temperatures. Although different 
aspects of these theories found good agreement with ex- 
periment in quantum dots 0] , a more general theory that 
accounts for the effects of temperature and dephasing at 
arbitrary fields and interaction strength remains to be 
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FIG. 1: (Left) Quantum dot with magnetic field B and dc bias 
voltages Vi,2 at the contacts and Vo at the gates with capaci- 
tance C; an additional lead (j> models dephasing. (Right) Nor- 
malized fluctuation of the (anti)symmetric component {Q a )Qs 
of the second order nonlinear conductance for a coherent dot 
as a function of flux $/<& c . Importantly, a crossover occurs at 
•frc <C <&o = eh/c, because of the long time an electron spends 
in the dot. 



developed. Such a theory is our main goal. 

This Letter presents a theory of fluctuations of conduc- 
tance nonlincarity in open chaotic dots within Random 
Matrix Theory (RMT). A key result of our work is il- 
lustrated in Fig. 2] The (anti) symmetric parts of the 
non-linear conductance strongly fluctuate from sample 
to sample due to quantum effects. These fluctuations are 
sensitive to the flux <!> through the dot. The fluctuations 
of the symmetric part Gs(&) decrease as the magnetic 
flux increases, while the anti-symmetric part Ga(&) has 
a stronger (linear Q) dependence at low fields. As the 
flux grows, the values of fluctuations reach their satura- 
tion values 0- The asymptotic values are equal for £? Q ,s 
as expected from the linear combination of uncorrelated 
contributions of random sign. 

The correct definition of the crossover scale $ c is im- 
portant for the quantitative comparison of the theory 
with experiment. It determines the slope of Q a at small 
flux $ (see Fig. QJ. One might naively expect <!> c ~ <f>o 
but importantly we find $ c -C $o in agreement with 
experiment Q . The scale of the crossover flux <I> C corre- 
sponds to a flux quantum <I>o = eh/c through a typical 
trajectory of an electron and not through the area of the 
dot. In a chaotic quantum dot the time Td an electron 
typically spends inside the dot is usually much larger then 
the ergodic time r erg necessary to explore its phase space. 
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During Td/r er g random attempts to explore the dot with 
flux <E> through its area, the flux penetrating the elec- 
tronic trajectory scales with $(r d /T org ) 1 / 2 . Therefore the 
crossover universally occurs at $ c ~ $o( T erg/' r d) 1/ ' 2 <C $ 
for ballistic or diffusive dots (the diffusive approach 
ignores Td 3> r crg , so $ c ~ $o of Ref. [j| corresponds to 
the flux quantum through the area of a diffusive dot). 

Below we find the fluctuations of Q a for arbitrary flux 
$, temperature T and capacitance C and compare in de- 
tail with previously considered limits 0, E} • We numer- 
ically investigate the effect of dephasing on CJ a at high 
magnetic fields, low temperatures, and strong interac- 
tion, which is relevant for experiments. Although the 
dephasing diminishes the conductance fluctuations 0,0, 
the uniform (locally weak) dephasing is found to have a 
stronger effect. We conclude with a comparison of theory 
and experiment. 

System. The 2D quantum dot, see Fig. ^ is bi- 
ased with dc voltages Vi,2 at contacts with N\jt ballistic 
channels, and by the voltage Vb at the gates with ca- 
pacitance C . This capacitance defines the strength of 
the Coulomb interaction in the dot (C — > corresponds 
to strong repulsion) |loj . The dot is in the universal 
regime when the Thouless energy i? T = ft/r crg is 
large, so that the results are applicable to dots with area 
A = irL 2 (taken circular), either diffusive with mean free 
path I <C L, or ballistic, with I 3> L and chaotic classical 
dynamics (in the latter case the substitution I — > ttL/A 
is used [l2j). The mean level spacing A = 2nh 2 /(m*A) 
and the total number of conducting channels N together 
define the dwell time Td = h/(NA) ^> r crg . A dephasing 
with rate 7^ = N ip TA/2tt is introduced with the dephas- 
ing probe model: a fictitious probe with N v channels of 
transparency T is attached to the dot 0. We also re- 
quire that eV -C NA and treat the nonlinearity only to 
{eV) 2 . Scattering is spin-independent and this spin de- 
generacy is accounted for by the coefficient v s . We use 
RMT for the energy-dependent scattering matrix S{e) 
and refer a reader to reviews flfl Hll| for details. 

The electric potential U in the dot is taken uniform. If 
the screening length is much larger than the Fermi wave- 
length, WKB can be applied. As a consequence, elec- 
trons with kinetic energy e have a well-defined electro- 
chemical potential e a — e — eV a in the contact a and 
i = e — eU in the dot. Therefore, transport depends on 
the Fermi-distributions f{e a ) and the scattering matrix 
S(e). The current in the contact a is I a = J del a {e) and 
for eV <C NA the spectral current I a {e) can be expanded 
in powers of eV: 



I a (e) = ^-^/(£ 5 )tr [l a S aS - l a $t(e)l a 5(e)] 
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In Eq. the total current I a is expressed in terms 
of the dimensionless linear conductance at energy e, 
9a@{ £ ) = tr(l Q £ Q/ 3 — l | g<S't ' (e)l a S(e)) and the nonlin- 
ear conductance g Q( 3 7 (e) related to d 2 I a /dVpdV 7 , which 
depends on U. To this accuracy U needs to be known 
only up to the first order derivatives, the character- 
istic potentials us — dU/dV$ [(J. The characteristic 
potentials us € (0, 1) are found self-consistently Q 
from current conservation and gauge-invariance require- 
ments and expressed in terms of the Wigner-Smith ma- 
trix Q = £^5/(2^) 0: 
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-J def(e)trQl s 
C/e 2 v s -f def'{e)tiQ 
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To leading order in N the mesoscopic average of Eq. 
is flux-insensitive. However, the fluctuations of us are 
strongly dependent on <!>, and determine the asymmetry 
of the non-linear conductance. These derivatives are used 
to express the conductances g a pj(s)'- 

3a/3 7 (e) = [8f3yg' a f)(£) ~ "/33a 7 0) ~ ^g' a p{e)]/2, (3) 

where the prime stands for energy derivative. The matrix 
S depends on magnetic field, 5($) = S T (— $), so that 
uo(— $) = uo(&), but importantly u$(&) lacks such sym- 
metry. Therefore quite generally g a p^($) 7^ 9ap-y{~^)- 
However some symmetries still hold for g a p^, which is 
revealed in the (anti)symmetric to $ — > — $ components 
of conductance {Q a )Q s (in units of inverse energy), 
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which we investigate now in detail. First we derive their 
dependence on temperature T, magnetic flux <E> and ca- 
pacitance C for coherent multi-channel dots, N 1, and 
later investigate partially coherent dots at T, C — > and 
high <I>. 

Coherent dot at arbitrary T, <E> and C. In a two- 
terminal dot without dephasing, we use gauge invariancc 
and set V2 = and consider derivatives with respect to 
Vi only. We define G a (s) = G a (s), ill and introduce a 
traceless matrix A = (N 2 /N)t 1 - (JVi/JV)l 2 such that 
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C/(e 2 v s )- J <fe/'(e)trQ 
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with fluctuating xi(e) = (A/27r)<9 e tr AS" 1 AS and x 2 (e) = 
(«A/27r)tr A[<S, d e S>]. The mesoscopically averaged G a 
vanishes, (Q a ) = 0, and we need to find correlations of 
Ga to leading order in N. To this end the products of 
S(E, $) and S'(E', $') are averaged, and the pair corre- 
lators, Cooperon Ce-E' and Diffuson T>e~e', are intro- 
duced as X £ = (N x - 2me/A)- 1 ,X = C,V, with the 
flux-dependent effective number of channels Nx 01 : 



N c 
N v 



= N ■ 



($± $') 2 hv F l 
JXK' 



(6) 
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The denominator of Eq. JBJ is a self-averaging quantity, 
((...) 2 ) = ((...)> 2 = (C/^A)) 2 , with the electrochemi- 
cal capacitance = C/(l + CA/(v s e 2 )) 17]. The func- 
tions xi{ £ >&) an d Xi{^ are uncorrelated, and their 
auto-correlations 
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readily allow one to find correlations of Q a ^ s at different 
magnetic fields and/or temperatures. In what follows we 
present the results for the variance of various experimen- 
tally measurable quantities. Calculated at the same T 
and <I>, they are given by very similar expressions. Most 
important are the mesoscopic fluctuations of S a ,s- 



T± = 
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The variance (Q 2 ) for a realistic dot |lj| with N\ — N 2 is 
given by Eq. (8) with T- — > J-+ . As a consequence, high 
magnetic field reduces the fluctuations of Q s by a factor 
two, see Fig. ^ The Eqs. (|8I9|) at high and low magnetic 
fields and weak interaction are of special interest and 
considered below in detail, since they allow us to compare 
with the results of Refs. [sEl- 

Low magnetic fields. When Nc ~ N the magnetic 
field is low and can not destroy time-reversal symmetry 
(TRS). In the limit T = 



A C 



NfN, 
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At small T < NA/2ir we have Q 2 a {T) / g 2 a {Q) w 1 - 
8(7r 2 T/A^A) 2 . For T > NA/2tt Eq. JEl is multiplied 
by 2(NA/24T) 2 and later used to compare with experi- 
ment . We point out that the TRS-breaking flux that 
destroys weak localization correction in open dots [lflj . 
obtained by other methods [lol Eoj l for chaotic and dis- 
ordered dots has the same dependence ^l] on as $ c 
in Eq. |TU|). 



For weak interaction, C/e A 



00, 



wc find 
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2irC h 2 $ r 
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We underline that Eq. Ijlljl holds for chaotic dots inde- 
pendently of the nature of scattering, diffusive or ballistic. 
The first term, rewritten via 2e 2 /C — ► /3/isA, reproduces 
the result of Ref. |jj if H/ta were substituted by the escape 
rate fi/r crg of an open diffusive sample (as is common for 
the crossover from ballistic to diffusive systems). 

The most important is the third term, large as 
Td/^crg ^> 1 for chaotic dots. It universally predicts the 
only relevant scale <f> c <C $0 rather then <E> C ~ $0 as 



(3 1 



> 



0.5 





- N d, =1 




N k =3 








- N =10 




... N =20 







0.5 



FIG. 2: var Q a , normalized by (71-/I6A) 2 , as a function of 
dephasing -y v for N = 2, Ni = 1 for iV^ = 1, 3, 10, 20. 



stated in Ref. ;S] . Data of Ref. 0, where a nonlinearity 
with $ sets in at flux $ <C 'I'o, substantiate our estimate. 

High magnetic fields. When N c > N the TRS is fully 
broken, C/V->0, and at T = the functions T± = 1/N 2 
reproduce the result of Ref. 8- At T > NA/2n the 
asymptotes are T±{T) » A/{\2TN). For T = and 
weak interaction we obtain 



ilziV 

ACh 2 ) 



TT 2 N 6 ' 



(12) 



where the first term reproduces the result of Ref. 
if H/rd — > h/r crg , and the second fully accounts for 
a possible asymmetry in the contacts. Therefore, at 
$ ~ $0 > & c the result coincides with that of Ref. 
up to a numerical coefficient. 

Partially coherent dot. Dephasing with rate 7^ is 
treated using the dephasing probe model : the current 
into the dephasing probe tp is zero, I v {e) = at every 
energy e. The probe generates currents in the leads a 
due to a "voltage" V v (e) at the probe. (V v (e) defines the 
distribution f(e — eV ip (ej) at the probe). For simplicity 
we take here a strongly-interacting dot, C = 0, and find 
v s (s) = dV tp (s)/dVg = -g v s{£)/9 w {£) and u s : 
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J def'(e)tr(Ks + -k v v 5 )Q 
Jdef'(e)trQ 



(13) 



For a dot at T = and $ 3> <I> C we numerically consider 
the antisymmetric component Q a : 



Qa 



(ffii + 9i v v+ + g' vl v- + g' v+v-) , (14) 



with u± = iti(±$),u± = ui(±$). If the dimensionless 
dephasing rate 7^, normalized by 2tt/A, is fixed, one 
can vary the transmission T together with the number of 
channels N v to go from uniform dephasing, N v ^> 1, to 
non-uniform dephasing with a small number of perfectly 
conducting channels in the probe |2lJ. We follow Refs. 
[22I |23| to generate Q and S for the broken TRS for a 
non-ideal coupling with the probe. The results for var Q a 
as a function of 7^ for N = 2, Ni = 1 are presented in Fig. 
El We expect that uniform dephasing damps Q a stronger, 
which is clearly seen in Fig. [21 (the uniform limit N v ^> 1 
is reached at N v ~ 10). In the limit 7^ ^ TV we have 
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G a — > 0. This behavior of £?(7p) could also be explored 
in experiments with a real additional probe (in the limit 
eV, T, A <§; NA the results of the dephasing probe |l3| 
and the inelastic probe [24| models coincide). 

Comparison with experiment. Zumbuhl et mea- 
sure the statistics of the (anti)symmetrized with respect 
to magnetic field B conductance {gB-)gB+ at various 
V, B. At low eV they correspond to (2v s e 3 /h)VG a 
and (v s e 2 /h){gn + 2eVG s )- Of particular interest are 
their rms 8gs± and the coefficient 6a — 8gs- /(VB) at 
eV, B — > 0. The measurements of da are performed at 
$ <C $o and T = 4/zeV for samples with A = TjieV 
for N = 2,4,8. Theory (see Q or Eq. © EJ) for 
$ = $o) predicts that the scaled coefficient 8a' (N) = 
SaN 2 Q hA/(2e 3 A) is independent of N. Experiment 
indeed finds a reasonable agreement with this prediction 
only at N = 8 where the measured value is 6a' (8) = 1.1 
versus the predicted 6a' th = TtC^/2C. Instead of a con- 
stant 6a' the experiment revealed a strong dependence 
of 6a' (2) re 0.015, <W(4) w 0.36 on iV. Our theory pro- 
poses an explanation of this unexpected growth. Inter- 
estingly, it turns out that, since TYd ~ h, we can use 
the low-field high-temperature asymptote of Eq. ©. It 
provides a behavior similar to the experiment: 8a' (N) = 
(ir 2 C^/24CT)(NAhv F /L) 1 / 2 . At Hv F /L ~ 250/ieV, one 
finds 6a' (N) pa 4.5Va"C m /C growing with AT but not 
nearly as steeply as observed. 

Experiment P| derives dephasing rate from weak lo- 
calization measurements, 7^ = 0.3. Numerics performed 
at <I> 3> $ c shows that uniform dephasing 7^ = 0.3 
diminishes 8a' (N) only by a factor 0.4, 0.65 and 0.85 
for N — 2, 4, 8 respectively. This is close to the result 
one would obtain with the semi-empirical substitution 
N — > N + j v in Eqs. JH1 EJ, which would strongly di- 
minish fluctuations in the few-channel dots. As a result, 
assuming non-ideal screening, C^/C < 1, and strong de- 
phasing one can fit our results to the experimental data 
of Ref. 0- If the experimental estimate of 7 y = 0.3 and 
the assumption that C\ t /C — 1 are used, our analytical 
results and experimental data disagree. Therefore both 
an accurate determination of the dephasing rate and an 
independent measurement of the capacitance are needed. 

We would like to point out that the symmetrized part 
of the non-linear conductance gs+ permits just such 
an independent determination of the capacitance ratio 
Cu./C. For symmetric (Aq — N2) coherent realistic dots 
[l8j at small eV it holds 

V = var deV 9B+ ( Vse 2 A \ 2 
\C) (var 5B+ )2 { h 2tt) ■ [ > 

at arbitrary flux <I> and temperature T. Since C^/C is 
apriori unknown, Eq. (|15f) presents an independent way 
to measure screening. Independent capacitance measure- 
ments and in particular also experiments for asymmetric 
dots {N\ ^ iVjj), or measurements of the nonequilibrium 
distribution |25j are needed for a more detailed compar- 



ison of theory and experiment. 

Conclusions. The Onsager-Casimir relations are a cor- 
nerstone of irreversible linear transport. We go beyond 
the linear regime and consider the quantum fluctuations 
of transport properties due to finite voltage and inter- 
actions by the example of open chaotic dots. Our key 
universal result, confirmed by experiment Q, is that the 
only relevant magnetic-flux scale is not the flux quantum 
through a dot $0 but rather <I> C <C $0 determined by 
a long dwell time of electrons inside the dot. Since our 
theory accounts for a wide range of parameters (tem- 
perature, flux, dephasing, contact widths), it provides a 
basis for further experimental investigation of the non- 
equilibrium transport. 
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